IEEE Transactions on Biomedical Engineering
● Institute of Electrical and Electronics Engineers (IEEE)
Preprints posted in the last 90 days, ranked by how well they match IEEE Transactions on Biomedical Engineering's content profile, based on 40 papers previously published here. The average preprint has a 0.04% match score for this journal, so anything above that is already an above-average fit.
Shenbagam, M.; Venkataraman, S.; Mukherjee, B.
Show abstract
Non-invasive human-machine interfaces (HMIs) are critical in developing prosthetic systems that offer intuitive, simultaneous, and proportional control over multiple degrees of freedom (DOFs). This study introduces a novel system for intuitive concurrent control of hand and wrist movements using sonomyography based imaging of muscle activity. Our method uses a sparse set of ultrasound scanlines to reduce computational complexity while enhancing usability. We evaluated four regression techniques for wrist and hand angle prediction, focusing on performance with a reduced sonomyographic feature set. We also explored the feasibility of a sonomyography-based system by simulating various factors that could affect prediction, including feature selection and scanline count. Our findings demonstrate that Gaussian process regression excels in predicting wrist and hand angles with just eight equispaced transducers in offline settings. Real-time evaluations with 10 non-disabled participants showed a 93 % success rate for two-DOF tasks using linear regression. The system was tested with an individual with amputation, achieving a 46 % success rate for two-DOF control in a 2D space, even though the ground truth data for model training was collected from the contralateral limb. This study validates our sonomyography-based approach for accurate wrist and hand angle estimation, reducing complexity and demonstrating potential in real-world scenarios.
Du, Z.; McManus, L.
Show abstract
Objective: Decomposition of high-density electromyographic signals enables non-invasive analysis of individual motor unit (MU) behavior, but reliable interpretation of physiological changes in health and disease depends on accurate MU discharge detection. This accuracy is compromised by mixed source estimates, where high amplitude peaks are associated with discharges from more than one MU. We introduce a post-decomposition framework to identify and separate suspected mixed sources using PCA-guided source refinement. Method: For each suspected mixed source, extended and whitened EMG vectors were extracted at source peaks and projected into a low-dimensional PCA subspace. This subspace highlighted MU-specific differences across candidate discharges, including subtle or spatially localized features of the spatiotemporal MUAP profile. Clusters in the PCA subspace were used to initialize source estimates for the constituent MUs. During iterative source refinement, source peak amplitudes were reweighted according to the distance of their corresponding points from the associated cluster center. Particle swarm optimization selected the reweighting factor that minimized the coefficient of variation of inter-spike intervals (CoVISI). Results: The algorithm separated mixed MU sources in simulated and experimental HDsEMG data. In simulated data, resolving mixed sources increased median rate of agreement (RoA) by >40%. In experimental recordings, MU yield increased by 1.27 per trial and CoVISI decreased by 0.28 (33% RoA improvement). Conclusions: PCA-based representation enhanced separability between MUs with similar MUAP profiles, while distance-based amplitude reweighting reduced re-merging during source refinement. Significance: This framework resolves merged MU discharge trains, improving decomposition accuracy and recovering MUs that might otherwise be excluded by quality thresholds.
Yuan, Y.; Li, W.; Zhu, L.; Su, H.; Yu, H.; Wang, H.; Lin, G. N.
Show abstract
Freezing of gait (FoG) in Parkinson's disease is a brief but hazardous gait failure that often precedes falls. For wearable cueing or other closed-loop assistance, a detector that reacts only after FoG onset is usually too late; the more useful task is to recognize the pre-freezing transition from physiological signals. This study presents PreFoGNet, a dual time-frequency deep learning framework for early FoG prediction using plantar pressure signals. The temporal stream combines a multi-scale Inception encoder with a bidirectional Mamba module to capture both short contact-related transients and several-second gait deterioration without the quadratic cost of attention. In parallel, the frequency stream uses band-wise spectral modeling and attention-based gating to emphasize physiologically meaningful changes in the locomotion, freeze-related, and high-frequency bands. On the WearGait-PD dataset, with a 2 s prediction horizon and subject-wise evaluation, PreFoGNet achieved a sensitivity of 93.94%, a specificity of 89.76%, a G-Mean of 0.9183, and an AUC-ROC of 0.9607. It outperformed classical machine-learning and deep learning baselines, and retained usable performance under moderate noise and single-channel loss. Additional horizon analysis showed that plantar pressure contains a stable pre-freezing signature within 0-3 s before onset, with a practical prediction boundary of approximately 6-7 s. These findings suggest that time-frequency modeling of plantar pressure is a promising signal-processing route for wearable FoG early-warning systems.
Hassan, M. W.; Crook, K.; Gi, Y. J.; Lee, J.; Hossain, M. M.
Show abstract
Objective: This study aims to develop and validate a quantitative, depth-resolved anisotropy imaging framework that extends ARFI-based focal degree-of-anisotropy (DoA) estimation into two-dimensional mapping by modeling the depth-dependent relationship between shear modulus ratio (SMR) and peak displacement ratio (PDR). Methods: We propose APRIL (Adaptive Polynomial Regression for anisotropy Imaging via ARFI-induced DispLacements), a framework for quantitative, depth-resolved DoA imaging that adaptively selects polynomial regression or shape-preserving spline interpolation based on excitation PSF asymmetry. Training data were generated using an LS-DYNA3D + Field II simulation pipeline in homogeneous transversely isotropic media (SMR 0.9-4.9). Testing included shifted SMRs under varied acoustic conditions and three heterogeneous inclusion configurations (anisotropic inclusion in isotropic background and vice versa). Experimental validation was performed in an in-vivo murine tumor model over the time, ex-vivo chicken breast, and tissue-mimicking gelatin phantoms, using a Verasonics system with an L11-5v transducer. Results: APRIL achieved depth-resolved SMR prediction errors below 9% over 10-30 mm, with highest accuracy in the focal region (MAE 2.3%, RMSE < 0.1) and stable performance across PSF transition zones. In heterogeneous phantoms, it reconstructed anisotropy maps with SSIM up to 86% and MPE below 7%, accurately delineating inclusion boundaries. Under acoustic parameter variations, mean absolute errors remained below 10%, demonstrating robustness to system and tissue heterogeneity. Conclusion: APRIL enables robust, two-dimensional anisotropy imaging beyond focal estimates. Significance: The method provides a physically grounded and generalizable framework for clinically viable anisotropy biomarkers in muscle, tendon, kidney, tumor and breast tissues.
Trisha, S. M.; Rahman, M. A.; Hassan, M. W.; Gi, Y. J.; Lee, J.; Hossain, M. M.
Show abstract
Viscoelastic characterization of tissue has significant diagnostic value in oncology, as tumor progression alters both elasticity and viscosity in ways that neither property alone can fully capture. Existing acoustic radiation force (ARF)-based methods such as Viscoelastic Response (VisR) ultrasound estimate relative elasticity and viscosity through per-A-line nonlinear model fitting, which is computationally intensive and requires auxiliary simulations to correct elasticity-dependent bias. This work presents VESTA (Machine Learning-Enabled Estimation of ViscoElastic Ratios from On-Axis Spatio-Temporal ARFI Features), a two-stage data-driven pipeline that predicts elasticity ratio (ER) and viscosity ratio (VR) directly from seven normalized ARFI displacement features at the A-line level, without model fitting or compensation. Stage~1 is an MLP classifier that detects inclusion boundaries from normalized peak displacement and negative peak velocity ratios; Stage~2 is a dilated Conv1D regression model that estimates ER and VR along the full axial sequence using the predicted mask alongside displacement features. The pipeline was trained on 500 simulated inclusion scenarios spanning three geometries, five focal depths, two F-numbers, and a broad range of material contrasts. In silico, mean predicted ER and VR were within 12\% of ground truth across all geometries, with performance best when ER and VR were moderate or decoupled. Experimental validation on a chicken breast phantom demonstrated plausible generalization to real tissue heterogeneity. Applied to an in vivo murine 4T1 breast cancer model, the pipeline tracked treatment-related attenuation of mechanical contrast in paclitaxel-treated tumors relative to controls over a 36-day imaging period, supporting its relevance for tumor monitoring.
Chen, B.; Subramanian, S.
Show abstract
BackgroundObjective tools for longitudinal dietary monitoring, despite its importance in the health triad of diet, sleep, and exercise, remain limited. Widely available ambulatory options include continuous glucose monitoring, which is a delayed response, and manual logging, which is rife with human error. Clinical tools such as gastric emptying scintigraphy are impractical for everyday use. High-resolution electrogastrography (HR-EGG) offers an alternative by treating gastric myoelectric activity as a biomarker of digestive state. However, its utility for ambulatory meal detection remains unclear. We hypothesize that HR-EGG and accelerometry together encode distinct postprandial gastric signatures to enable automated meal detection, and that postural context represents a relevant source of variation in signal detectability. MethodsHR-EGG and accelerometry data were collected from seven healthy adults across sixteen 150-minute meal sessions under IRB-approved protocol. Each session included a 30-minute fasted baseline, consumption of a standardized meal, and 90-minute postprandial period of sitting, walking, and lying in a randomized order. Features of the gastric slow wave, including raw and normalized bandpower, phase gradient directionality (PGD), wave direction, and wave speed, were extracted alongside triaxial accelerometer magnitude. A dilated one-dimensional convolutional network (1D CNN) was trained to classify meal consumption at five-minute resolution using leave-one-subject-out cross-validation. Postural effects on gastric myoelectric metrics were assessed using the Friedman test. ResultsThe model achieved a mean AUROC of 0.925 (95% CI: [0.857, 0.993]) and mean AUPRC of 0.824 (95% CI: [0.668, 0.980]; null model: 0.20). Feature ablation showed PGD as the most informative input ({Delta}AUPRC = -0.188), with wave propagation speed the least informative ({Delta}AUPRC = -0.105). Walking produced the highest signal-to-noise ratio (9.94 dB), lying had the most stable gastric rhythm (89.1% normogastric), and sitting demonstrated the greatest frequency instability (dominant frequency standard deviation = 0.825 cpm). ConclusionA dilated 1D CNN applied to spatiotemporal HR-EGG features enables temporally aware passive meal detection across ambulatory contexts. This study framework addresses a gap between clinical need for objective dietary monitoring and the limitations of current detection methods.
Ghaffarzadeh, P.; Chakraborty, D.; Aslansefat, K.; Dostan, A.; Papadopoulos, Y.
Show abstract
Ground reaction force (GRF) measurement remains largely confined to instrumented laboratories, limiting longitudinal monitoring in daily life. This article presents an edge-first wearable system for estimating vertical GRF from consumer smartwatches. Two Apple Watch Series 6 devices worn at the wrist and waist stream 12-channel inertial data at 100 Hz to an iPhone, where preprocessing, storage, and inference occur locally without cloud dependence. The proposed GRFNet-MultiScale model is a compact temporal convolutional network with four dilated residual blocks and a global context branch. Under leave-one-subject-out evaluation on 539 stance windows from 10 healthy participants, the dual-sensor system achieved a mean Pearson correlation of 0.798 with an RMSE of 257 N, while a wrist-only configuration retained 82.5% of dual-sensor correlation. Temporal attribution remained stable across validation folds and identified early-stance wrist acceleration as the dominant reproducible signal. The system is strongest for cyclic locomotion.
Das, S.; Sharma, K.; Sarkar, S.; Gonsalves, K.; Srinivasan, U. S.; VARMA, H.
Show abstract
SignificanceDiffuse Correlation Spectroscopy (DCS) is an established technique for non-invasive monitoring of Cerebral Blood Flow (CBF), but existing applications primarily measure CBF changes averaged over cortical tissue volumes. A method capable of targeting blood flow within a particular intracranial artery would expand the utility of DCS for vessel-specific cerebral perfusion monitoring and continuous bedside assessment. AimWe aim to investigate the feasibility of targeted, non-invasive monitoring of blood flow in the Anterior Cerebral Artery (ACA) using DCS through optimization of probe geometry and placement. ApproachA custom-built DCS system operating at 785 nm was used to probe ACA from the glabellar region. Source-detector (SD) separation, probe orientation, and probe location were systematically optimized using lower-limb motor tasks and a mental arithmetic task. The optimized configuration was evaluated using an ACA-mimicking multilayer phantom and validated in forty healthy volunteers during lower-limb activation tasks and postural changes. ResultsAn SD separation of 17 mm, vertical probe orientation, and glabellar placement provided the highest sensitivity to ACA-related blood flow changes. Phantom experiments demonstrated sensitivity to flow changes in a vessel located 45 mm beneath the scalp and showed that the frontal sinus, cerebrospinal fluid, and the absence of cortical tissue beneath the glabella along the longitudinal fissure together provide the best optical window for probing deep ACA flow. Using the optimized probe configuration, significant increases in relative CBF were observed during standing leg marching (66.4 {+/-} 38.6%) and supine leg crunches (39.4 {+/-} 32.2%) (p < 0.01), with consistent responses during supine-to-stand postural transitions. ConclusionThe proposed DCS approach enables targeted measurement of blood flow within the ACA territory. This technique provides a framework for continuous, non-invasive monitoring of ACA perfusion and has potential applications in cerebrovascular monitoring and stroke care.
Liu, D.; Dutta, A.; Nadig, S.
Show abstract
The features of the PPG (photoplethysmography) morphology are known to reflect age-related cardiac and vascular changes. In most contemporary wearables, PPG signals are acquired from distal sites such as the wrist and finger. The superficial temporal artery (STA), accessible at the temple region, is reached via a shorter arterial path from the aortic root than the radial circulation, and may therefore carry hemodynamic and aging information with less distance-dependent attenuation. We hypothesized that the morphology of the PPG at temple region (STA) would show stronger and more numerous age correlates than the PPG at the wrist. To test this, we extracted a common set of 89 pulse-morphology features, spanning raw-waveform timing/amplitude/area measures, ratios among them, derivative-based ratios, and spectral harmonic-ratio features. We compared an in-house temple-worn device which has PPG as one of the sensors, with a publicly available Microsoft Aurora-BP wrist-worn PPG dataset, and tested each feature's association with age. We identified 14 robust age correlates at the temple region, compared to 3 at the wrist. The temple's correlates spanned multiple morphological categories and showed a larger age-association than at the wrist. These results support the hypothesis that the temple region may be a more robust PPG measurement site than the wrist to extract age-related cardiovascular information, which motivates further investigation of temple-based cardiovascular sensing.
Tecchio, P.; Schlaffke, L.; Bolsterlee, B.; Hahn, D.; Raiteri, B. J.
Show abstract
Muscle architecture shapes muscle function and changes with age, growth, training and disease, yet quantifying three-dimensional (3D) muscle architecture in vivo remains challenging. We introduce a hybrid fascicle tractography approach for freehand 3D ultrasound data that accurately reconstructs 3D muscle fascicles with respect to an objective, anatomically relevant coordinate system defined by the muscle's central aponeurosis. The hybrid approach combines Hessian-based fascicle detection with wavelet-based refinement to generate volumetric fascicle orientations. In a synthetic dataset with known ground truth, fascicle orientations and lengths were estimated with errors of [≤]2{degrees} and ~1.5%, respectively. In vivo, the approach detected physiologically plausible fascicle lengthening in the human tibialis anterior following a passive plantar flexion rotation, whereas diffusion tensor imaging of the same muscle did not. The proposed method enables anatomically relevant, objective and non-invasive quantification of 3D muscle architecture in vivo, providing a practical framework for applications in clinical and applied muscle physiology.
Lim, J.; Islam, R.; Raghavan, D.; Omofojoye, B.; Rodriguez, A. D.; Kiarashi, Y.; Hershenberg, R.; Clifford, G. D.; Kwon, H.
Show abstract
Mild cognitive impairment (MCI) is a clinically important stage preceding Alzheimer's disease and related dementias, in which cognitive and balance functions are commonly evaluated using standard clinical assessments such as the Montreal Cognitive Assessment (MoCA) and Mini-Balance Evaluation Systems Test (Mini-BESTest). These assessments are administered episodically by clinicians and may miss functional changes during everyday movement. Recent studies and prior work in the Charlie and Harriet Shaffer Cognitive Empowerment Program (CEP), a therapeutic environment supporting lifestyle intervention and naturalistic social interaction, suggest that wearable and passive behavioral sensing can monitor movement patterns associated with cognitive and balance function in older adults with MCI. However, it remains unclear whether passive waist-mounted IMU data collected during naturalistic movement and social interaction can quantify clinician-rated cognitive and balance outcomes, particularly at the subdomain level, in an interpretable and demographically fair manner. To address this gap, we analyzed weekly IMU recordings collected over 6 months from 44 older adults with MCI in the CEP and trained tree-based ensemble regression models to estimate MoCA and Mini-BESTest total and subdomain scores, with interpretability and demographic fairness evaluation. Our models achieved RMSEs of 3.677 for MoCA and 3.672 for Mini-BESTest, benchmarked against Minimal Detectable Change and Minimal Clinically Important Difference thresholds. Feature importance analysis showed distinct movement signal properties across assessments, with general movement intensity features most informative for MoCA and temporal gait features led by cadence most informative for Mini-BESTest. Demographic bias analysis identified sex-related model bias, mitigated through post-processing while maintaining performance. This study supports the feasibility of wearable-based estimation of clinical assessment scores in older adults with MCI during naturalistic activity, with comparable performance between sexes after bias mitigation. This advances the validation of passive sensing for home monitoring to support clinical decision-making and personalized interventions.
McCorkendale, B.; Rodriguez, R.; Fink, R.; Moore, M.; Romero, S.; Esmailie, F.
Show abstract
PurposeMild therapeutic hypothermia (MTH) preserves cochlear function in animal models and is now entering early-phase human trials for hearing preservation. However, the extent to which the human cochlea can actually be cooled, and the mechanisms underlying MTH, remain unclear, in part because blood perfusion is expected to oppose localized cooling. In this study we evaluated the impact of blood flow on human cochlear temperature exposed to the MTH device using a combined experimental and computational approach. MethodsTemperature measurements were obtained from a human cadaver skull exposed to a commercial MTH device. These data were used to validate a three-dimensional bioheat transfer model incorporating realistic skull anatomy. The validated model was subsequently extended to include physiological blood perfusion in the internal carotid artery; a major heat source located near the cochlea. Finally, the in silico model was further expanded to incorporate the surrounding skin and brain tissues. ResultsIncorporating blood flow in internal carotid artery substantially altered predicted cochlear temperature distributions, highlighting the importance of localized vascular heat transport in the human cochlea during MTH. Although cochlear cooling was attenuated in the presence of perfusion, the therapeutic effects of MTH may not depend solely on the magnitude of local intracochlear temperature reduction. Additional mechanisms, such as reduced facial surface temperature, may also contribute to its efficacy. ConclusionThe validated in silico model provides a physiologically realistic framework for evaluating human cochlear thermal responses, investigating MTH mechanisms, and optimizing temperature-based strategies for hearing preservation.
Alipour, A.; Acikel, V.; Gokyar, S.; Algin, O.; Oto, C.; Balchandani, P.; Demir, H. V.; Atalar, E.
Show abstract
PurposeTo enhance the SNR in MRI within a localized region of interest using a novel interventional wireless RF resonator probe combined with a dual-drive pTx system. MethodsA dual-drive body birdcage coil was operated in a linearly-polarized mode to decouple a passive RF resonator probe from the transmit field while maintaining the resonator in a receive-only coupled mode. The resonator was fabricated using standard microfabrication techniques and tuned to the Larmor frequency of a 3T MRI system. The 10-g specific absorption rate (SAR) distribution was simulated to identify potential hot spots around the resonator prior to heating experiments. To evaluate the interaction between the resonator probe and the linearly-polarized transmit field, SNR and flip-angle distributions were measured in a phantom. In vivo imaging studies were subsequently performed using the resonator probe in conjunction with the linearly-polarized dual-drive birdcage coil. ResultsTemperature measurements demonstrated a normalized temperature increase of less than 0.10{degrees}C, corresponding to a SAR value below 1.21 W/kg. Experimental flip-angle mapping confirmed effective magnetic decoupling of the resonator probe using linearly-polarized dual-drive transmission. An SNR enhancement factor of 1.6 was achieved within the region of interest in phantom experiments. In vivo imaging demonstrated a 2.0 {+/-} 0.2-fold SNR enhancement in the vicinity of the resonator probe. ConclusionA novel interventional approach for localized SNR enhancement in MRI was demonstrated using a wireless RF resonator probe and a dual-drive pTx system. The proposed technique enables local signal enhancement while minimizing transmit-field interactions, thereby facilitating safe interventional MRI and potentially improving image quality and diagnostic performance.
Fu, J.; Zhang, S.; Huang, H. J.; Rakhshan, M.; Wen, Y.
Show abstract
Motor unit (MU) decomposition using high-density surface electromyography (HD-sEMG) has been widely used to characterize MU behavior in neurophysiology and to build neural-machine interfaces for wearable robots. Recently, many open-source software tools for MU decomposition have been made available on GitHub, which could reduce the effort of researchers in the field. However, the consistency among these open-source tools has never been studied, making researchers hesitate to use them. In this study, we collected 7 open-source software tools on GitHub and applied them to decompose MUs from an open-source HD-sEMG dataset (including 11 isometric contraction trials) to investigate the consistency among these tools. To create a comprehensive MU pool for reference, we combined all unique MUs identified by seven tools, visually inspected and removed bad MUs, and manually edited all remaining MU spike trains. Across 7 tools for 11 trials, the number of identified MUs ranges from 167 to 736. The number of valid MUs after expert inspection ranges from 29 to 210, which is 10% to 72% of the reference pool. The rate of agreement between the raw MUSTs and the manually edited MUSTs ranges from 0.86 to 0.94, and the averaged number of edits per MU to correct misalignments ranges from 14 to 39. The results show inconsistency in the implementation and procedures of each tool, which results in an inconsistent number of identified MUs and valid MUs (29 vs 210). In general, a substantial amount of effort is required to process the raw MUSTs from each tool to conduct further research analysis. This study provided a guideline for using open-source software tools for MU decomposition and indicated that it would be beneficial to develop tools to automatically edit the MUSTs.
Qiu, C.; Li, D.; Huo, H.; Mishra, A.; Li, C.; Yin, K.; Wang, N.; Chen, J.; Yao, R.; Margolin, E. J.; Lipkin, M. E.; Zhong, P.; Ni, X.; Yao, J.
Show abstract
Urinary stone disease is a common urological condition with increasing incidence, particularly in developed countries. Laser lithotripsy (LL) has become a preferred minimally invasive treatment due to its high precision and low tissue damage. Recent studies suggest that cavitation plays a critical role in stone damage during LL, and three-dimensional passive cavitation mapping (3D-PCM) has emerged as a promising tool for detecting these events. However, clinical translation of 3D-PCM remains challenging due to limitations in imaging depth, field of view (FOV), and procedural compatibility. Here, we present a large-FOV dual-modality imaging system (3D-PCM and B-mode ultrasound) based on a large-aperture planar ultrasound array. Through array optimization and model-based reconstruction, our system achieves an expanded FOV of ~40*40mm^2 at a clinically relevant imaging depth of ~110mm, while maintaining high spatial resolution of ~0.6 mm laterally and ~0.4 mm axially. In vivo experiments in a porcine model demonstrate that the reconstructed cavitation distribution correlates well with stone damage. Our technology has the potential to provide real-time treatment feedback during LL without disrupting the standard workflow.
Ilovitsh, T.; Shapiro, G.; Gershman, Y.; Bismuth, M.
Show abstract
This study presents the use of sub-micron nanobubbles (NBs) as contrast agents for ultrasound localization microscopy (ULM), a super-resolution imaging technique that visualizes microvascular structure and flow beyond the acoustic diffraction limit. While ULM has traditionally relied on micron-sized microbubbles (MBs), the reduced dimensions and prolonged circulation times of NBs make them attractive candidates for localization-based imaging. However, their weaker acoustic responses present significant challenges for reliable detection and tracking. To address this challenge, we developed the ULM Master GUI, an interactive framework for optimization of the complete ULM processing pipeline. Using custom ultrasound-compatible wall-less gelatin flow phantoms containing vessel-mimicking channels and bifurcations ranging from 100 to 500 m, we demonstrate that NB-based ULM achieves velocity reconstruction and flow partitioning measurements comparable to conventional MB-based ULM. Across all investigated geometries, NBs faithfully reproduced the underlying flow patterns and hemodynamic behavior despite their substantially reduced acoustic scattering. These findings establish the feasibility of NB-based ULM, expand the range of contrast agents available for localization microscopy, and provide a foundation for future super-resolution ultrasound imaging using nanoscale acoustic contrast agents. The ULM processing GUI is publicly available at https://github.com/grisha1998/ulm-super-resolution-toolbox.
Mlynczak, M.; Rosol, M.; Korzeniewski, K.; Gasior, J. S.
Show abstract
Background and ObjectiveAccurately parameterizing dynamic, time-varying interactions in physiological systems is a methodological challenge, as global causal discovery methods may obscure transient, local fluctuations. This study introduces tempord, an open-source Python library designed to estimate local temporal orders and evaluate the short-term stability, directionality, and strength of causal links in non-stationary biological signals. MethodsThe algorithm estimates temporal relationships by keeping one signal stationary while iteratively shifting another one within a sliding window. To parameterize optimal inter-signal shifts (causal vector, CV), the framework utilizes linear modeling or time series distance metrics. The methodology was validated through a simulation study on synthetic bivariate signals with mathematically imposed dynamic phase delays, under both deterministic and noisy conditions. Furthermore, in-vivo capabilities were demonstrated by evaluating cardiorespiratory coupling dynamics across spontaneous and music-induced relaxation breathing states. ResultsThe simulation study demonstrated that the extracted CV trajectories precisely aligned with ground-truth temporal delays, assessed using mean absolute error and root mean square error for both noise-free and noisy synthetic data. In-vivo application demonstrated dynamic temporal stability and the detection of minor step changes during autonomic nervous system state transitions. ConclusionsThe tempord Python package bridges the gap between global causal discovery and local beat-by-beat statistical parameterization. It provides a robust "bottom-up" analytical instrument for investigating the transient mechanisms governing complex biological networks.
Lo, H. U.; Gao, Z.; Loi, H. F.; Cheng, S. K.
Show abstract
Surface electromyography (sEMG) is the most practical non-invasive interface for myoelectric prostheses, exoskeletons, and rehabilitation systems, but power-line interference (PLI) contamination and excessive digital pipeline group delay still limit its clinical adoption. This paper proposes a co-designed analog-digital correction system combining a high-CMRR front-end with an exponentially-windowed RMS (EMRMS) envelope estimator and a recursive single-tone PLI canceller. We present a closed-form CMRR model capturing the electrode-skin imbalance, and provide a complete stability analysis of the LMS canceller. The EMRMS estimator reduces the computational overhead from[O] (L) to strictly[O] (1) in both time and space complexities. Featuring no data-dependent branching, the algorithm achieves deterministic algorithmic execution time (zero jitter under an RTOS environment) and is natively compatible with fixed-point arithmetic on microcontrollers lacking a hardware Floating-Point Unit (FPU). A reference implementation reaches an 8.2 {micro}s median per-sample latency, yielding an end-to-end delay of[~] 30 ms -- leaving a generous >90 ms budget for electromechanical actuation -- while requiring an active CPU duty cycle of merely 1.6%, enabling prolonged deep-sleep intervals. Validation on the public Ninapro DB2 dataset demonstrates a 13.9 dB mean SNR improvement (averaged across 12 channels; single-channel comparison: 9.7 dB, Table 3) and a 70.0 {micro}V envelope RMSE against a length-200 rectangular reference. Paired Wilcoxon signed-rank tests confirm statistical significance (p < 0.001) over static baselines, and Pearson correlation analysis ({rho} = 0.993 {+/-} 0.0002) confirms strict morphological fidelity. The full open-source codebase and benchmarks are publicly released. O_TBL View this table: org.highwire.dtl.DTLVardef@299dc5org.highwire.dtl.DTLVardef@3519a0org.highwire.dtl.DTLVardef@2586aborg.highwire.dtl.DTLVardef@1ac5610org.highwire.dtl.DTLVardef@1465c46_HPS_FORMAT_FIGEXP M_TBL O_FLOATNOTable 3:C_FLOATNO O_TABLECAPTIONQuantitative comparison on a common 60 s segment of Ninapro-like synthetic sEMG (single channel) with a 3 mV 50.3 Hz mains tone slightly drifted from the static notchs design centre at 50.0 Hz, stress-testing the adaptive corrector under a frequency mismatch. The Ninapro multi-channel aggregate (13.9 dB) reported in Section 3.4 uses mains exactly at 50 Hz (matched notch) and so achieves a higher {Delta} SNR. "MAC/sample" excludes the EMRMS square root and the pre-computed LMS sine/cosine. C_TABLECAPTION C_TBL
Liu, Z.; Duncombe, P.; Napadow, V.; Handsfield, G. G.
Show abstract
Physiological cross-sectional area (PCSA) is defined as the summed cross-sectional area of all muscle fibers contracting in parallel and at optimal length. PCSA is widely used, yet the conventional equations used to compute PCSA were developed from two-dimensional (2D) interpretations of muscle architecture and may not accurately represent muscle fiber cross-sections in the case of real three-dimensional (3D) geometries of muscles. Related measures of functional cross-sectional area (FCSA) and geometric cross-sectional area (GCSA) were also developed and interpreted with simplified 2D representations. Using realistic 3D muscle architectures derived from medical imaging, we sought to investigate whether conventional definitions of PCSA, FCSA, and GCSA represent the summed cross-sectional areas of all parallel muscle fibers, the fundamental definition of PCSA. We found that none of these measures consistently represented this definition. Thus, we introduce the physiological cross-sectional surface (PCSS), a curved surface within a muscle volume that is everywhere perpendicular to the local fiber direction. We estimated PCSS in 3D muscle surface meshes reconstructed from MRI data, using fiber orientations derived from Laplacian fiber reconstruction. PCSS-derived estimates were compared with PCSA, FCSA, and GCSA across six muscles representing five architectural classes. PCSS differed from all conventional measures, with the magnitude and direction of disagreement depending on muscle architecture. PCSS-to-PCSA ratios ranged from 0.771 to 1.399, while GCSA underestimated PCSS by up to a factor of 2.256 in bi- and multipennate muscles and overestimated it in muscles with more uniform fiber arrangements. PCSS demonstrated high geometric fidelity (perpendicularity>0.987) and robustness to fiber density across a tenfold range (coefficients of variation 0.39-3.95%). These findings indicate that conventional cross-sectional area measures do not consistently account for all fiber cross-sections in parallel within realistic 3D muscle geometries. PCSS provides a geometrically rigorous alternative that may improve estimation of functional muscle capacity from subject-specific imaging data.
Gonzalez Nunez, J. G.; Sabri, S.; Kebria, P.; Crook, J.; Brattain, L.
Show abstract
This paper presents a two-stage pipeline for implicit feature engineering in time series-based physiological stress detection using electrodermal activity (EDA) signals. In the first stage, we forecast three descriptive statistics of future EDA signals over short horizons (3, 5, and 10 seconds) based on a 60-second context window. In the second stage, a lightweight linear classifier detects stress from these predicted statistics. We evaluate three forecasting architectures spanning the domain expertise spectrum: a domain-specific bidirectional long short-term memory (BiLSTM) recurrent neural network, zero-shot and fine-tuned variants of Amazon Chronos T5 time series foundation model, and the Tabular Prior-data Fitted Network (TabPFN) applied to engineered physiological features. Experiments on the publicly available Wearable Stress and Affect Detection (WESAD) dataset, comprising chest-worn multimodal physiological signals from 15 subjects under baseline and stress conditions, demonstrate that the domain-specific BiLSTM achieves the highest classification performance, with area under the receiver operating characteristic curve (AUC) values ranging from 0.913 to 0.962. TabPFN follows with AUC values of 0.853-0.869, while Chronos variants yield 0.528-0.744. Notably, models using predicted features consistently outperform those using oracle features derived from the true future signals--the theoretical upper bound--suggesting effective noise filtering through learned sequence representations. Chronos models quickly reach performance saturation regardless of training depth, highlighting challenges in tokenizing continuous physiological time series. The proposed approach advances implicit feature engineering for wearable stress monitoring by leveraging forecasting as a powerful inductive bias, thereby improving robustness and providing insights into the limitations of the foundation model for physiological signals.